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Abstract 

Nourdin et al. [9| established the following universality result: if a sequence of off-diagonal homoge¬ 
neous polynomial forms in i.i.d. standard normal random variables converges in distribution to a normal, 
then the convergence also holds if one replaces these i.i.d. standard normal random variables in the poly¬ 
nomial forms by any independent standardized random variables with uniformly bounded third absolute 
moment. The result, which was stated for polynomial forms with a finite number of terms, can be ex¬ 
tended to allow an infinite number of terms in the polynomial forms. Based on a contraction criterion 
derived from this extended universality result, we prove a central limit theorem for a strongly dependent 
nonlinear processes, whose memory parameter lies at the boundary between short and long memory. 


1 Introduction 

In Nourdin et al. Q, a universality result was established for the following off-diagonal homogeneous poly¬ 
nomial form 

Qk(N n , f n , X) := Y, f n (ii,...,ik)X il ...X ik , (1) 

l<ii,...,i k <N n 

where f n is a sequence of symmetric functions on Z^ vanishing on the diagonals ( f n {h ,..., ik) = 0 if i p = i q 
for some p y q), and X = (X \, X'i -...) is a sequence of standardized independent random variables, and N n 
is a finite sequence such that N n —> oo as n —> oo. 

The universality result says that if Z = ( Z\. Z-j....) is an i.i.d. standard normal sequence and Qk{N n , f n , Z) 
converges weakly to a normal distribution as n —> oo, then the same weak convergence to normal holds if 
Z is replaced by X, where X is any standardized independent sequence with some uniform higher moment 
bound. 

It is natural to try to eliminate the finiteness of N n in the preceding result. This extension was mentioned 
in Remark 1.13 of Nourdin et al. @, but was not explicitly done. One would encounter a number of difficulties 
if one were to extend the method of proof used for finite N n to N n = oo. We will note, however, that this 
extension can be easily achieved using a simple approximation argument. We find it valuable to have such 
an extension and the corresponding contraction criterion ('Theorem 12.611 since it can be directly applied to 
limit theorems in the context of long memory. 

We consider such an application in Section [3] where we suppose that 
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and where the function a(-) behaves essentially like a homogeneous function with exponent a. The resulting 
polynomial form <5fe(/iv) is then the partial sum of a stationary process. The exponent a is chosen in such a 
way that the corresponding stationary process lives on the boundary between short and long memory. We use 
the contraction criterion to prove that a central limit theorem holds but with the nonstandard normalization 
y/N In N. This delicate case seems difficult to treat otherwise. 

The paper is organized as follows. In Section [2j we state the and prove the extension of the universality 
result (Theorem 12.1 1) . and as a byproduct, a criterion for asymptotic normality (Theorem 12.61) . In Section 
EU we state the critical limit theorem obtained by applying the criterion. In Section 13,31 and 13.41 we give 
the proofs. 

2 Universality of homogeneous polynomial forms 

Let £ 2 (Z fe ), k > 1, denote the space of symmetric square summable functions on Z fc vanishing on the 
diagonals equipped with the discrete L 2 norm. Let X = (Xi, X 2 ,...) be a sequence of independent random 
variables satisfying E Xi = 0 and EX 2 = 1. By modifying the notation Jl]), one defines for / £ £ 2 (Z k ): 

Qk(f, X) := Y, 

— 00<ll ,...,2fc<00 


One has 

EQ fc (/,X) = 0. 

Consider now two homogeneous polynomial forms Qfci(/i,X)) and Qk- 2 (f 2: X), where f\ £ £ 2 (Z fel ) and 
f 2 £ £ 2 (Z k2 ). Then the covariance of and Qfc 2 (/ 2 ,X) is 

(A,/2> := EQfc(/ 1 ,X)Q fc (/ 2 ,X) (2) 

= f fc! E-oo<i 1 ,... > i k<0 c/i(U---,4)/2(U---,4), iik 1 = k 2 = k-, 

|o if h ± k 2 . 1 


We then have the following extension of Nourdin et al. Q Theorem 1.2: 

Theorem 2.1. For each j = 1 suppose that kj > 2, and let he a sequence of functions in 

£ 2 (Z fc i). Let E be an m x m symmetric non-negative definite matrix whose each diagonal entry is positive. 
Assume in addition that 

sup E < 00. (4) 

n .. . 

— 00<Zi j <oo 

Then the following two statements are equivalent: 

1. For every sequence X = (Xi, X 2 ,...) where Xi,X 2 ,... are independent random variables satisfying 
EX; = 0,EX 2 = 1, and 

supE|Xj| 3 < oo, (5) 

i 

the following joint weak convergence to a multivariate normal distribution holds: 

(Q fc ,(/„, j; X))J =i 4x(0, E). (6) 

2. For a sequence Z = (Z\, Z 2 ,...) of i.i.d. standard normal random variables, the following joint weak 
convergence to a multivariate normal distribution holds: 

(Q fc3 .(/„,„Z))J =1 4X(0,E). (7) 


2 



Remark 2.2. Condition © can be re-expressed as 


supEQ fc . (/„j,Z) 2 = /cjisup ^2 < oo. ( 8 ) 

n n 

— <oo 

Remark 2.3. One can recover Nourdin et al. @] Theorem 1.2 from Theorem l2.1l bv replacing f n j(i i,... 
with fn,j (fl 7 • • • ) ^kj ) ll<ii <N n (il} • • • 3 ifcj ) • 

Remark 2.4. In the one dimensional case: m = 1, one can relax the assumption © by sup^ E|Xj| 2+<5 < oo 
for any 5 > 0. See Theorem 1.10 of Nourdin et al. [9j. 

Proof of Theorem \2.1\ We need to prove that © implies ©. Define the lV n -truncated functions 

fn,j (f 1? • • • ) fn,j\i Is • * • 5 ^kj ) 1{ — N n <ii<N n <iV n } > j = 1 ,,171. 

For any n £ Z+, we can hnd 7V„ large enough, so that for all j = 1,..., m, 


E 


Qkj ( fn,j , Z) - Qkj (fn,j , Z) 


= E 


Qkj (fn,j,X) - Q fej (/ njJ ,X )| 2 = kj\\\f n j - fnj< i. 


(9) 


Assume without loss of generality that N n —► oo as n —> oo. By © and ©, one has 

(Q kj (fnj, Z))™ i 4 X( 0 ,E). 

Using the original version of the universality result in Nourdin et al. !9J Theorem 1.2, one gets 

(Qfci(/n,i,x))™ i 4iv(o ) E). 

The conclusion © follows from © and m- 


( 10 ) 


□ 


Remark 2.5. Using the same argument as in the preceding proof, one can eliminate the finiteness of N n 
in © in the following related universality results for homogeneous polynomial forms: (a) Theorem 1.12 of 
Nourdin et al. [3] concerning for convergence to a % 2 distribution; (b) Theorem 3.4 of Peccati and Zheng 111 
which is the counterpart of Theorem 12.11 here with Zf s being standardized Poisson random variables. 

Theorem m gives rise to a practical criterion for the convergence ©. We first introduce the discrete 
contraction operator: for / £ £ P (Z P ) and g £ t q (fL q ), p, q > 2, we define 


OO 

(/ *r S)(* 1, • • • ! ip+q—2r) = 'y ) /O'l, • • • ) jri *1) • • • > ip—r^dijli ■ • • ) jr, ip—r+ 1, ■ ■ • i ip+q—2r ) (H) 

3 ii--->Jr=—oo 


for r = 0,... ,p A q, where in the case r = 0 it is understood as the tensor product. 

Theorem 2.6. Let {/n,j(')> n £ Z + } a sequence of functions in £ 2 {Z kj ) satisfying ©), j = 1 
where kj > 2. let E be an m x m symmetric non-negative definite matrix whose each diagonal entry is 
positive, such that 

£(bi) = lim (f n ,i,fn,j), (12) 

n— >oo 

where (■, •) is defined in m■ Then the following are equivalent: 

1. For every X = (Xi, X 2 ,...) with Xj’s being independent random variables satisfying E Xj = 0, EX 2 = 1 
and supjE|Xj | 3 < 00 , we have the following joint weak convergence to normal: 

(Q kj (/„,,-, X))J =1 4 JV( 0 ,E). (13) 
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2. The following contractions are vanishing: 


lim II f n j * r fn,jhki-2r = 0, for all r = 1,..., kj - 1 and all j = 1,..., to. 

n—>• oo J 


(14) 


where || • \\k denotes the discrete L 2 norm on £ 2 (Z k ). 

Proof. By Theorem l2.ll the statement 1 is equivalent to (Q kj ( f n ,j , Z))™^ —>■ N( 0, E), where Z is a sequence 
of i.i.d. standard Gaussian variables. Note also that each Q kj (fnj, Z) can be expressed as a fcj-tuple Wiener- 
Ito integral with respect to Brownian motion. For Wiener-Ito integrals, joint convergence to the normal is 
equivalent to marginal convergence, and marginal convergence is equivalent to the contraction relations. 
More precisely, by applying Theorem 6.2.3 and 5.2.7 of Nourdin and Peccati fl, one gets the equivalence to 
m- See also Theorem 7.5 of Nourdin et al. ||. □ 

Remark 2.7. We shall use the implication “Statement 2 => Statement 1” of the preceding theorem in 
the sequel. As for the reversed implication, namely, “Statement 1 => Statement 2”, the stipulation “For 
every” is important here, as well as in Theorem 12.II because there are random variables Xfs, for example 
Rademacher, that is X t = ±1 with probability 1/2 each, for which one may have convergence in (fl3l) even 
when (fl4l) does not hold (see Nourdin et al. [9|, Section 1.6, p.1956). 

Remark 2.8. One may wonder if the universality result extends to a continuous setting, namely, when 
Qk(fn) is replaced by a multiple integral on a Borel measure space (A,A,n): 


h(fn,0 


/ A k 


f n (x i,..., x k )^(dx i)... f(dx k ), 


where / € L 2 (A k ), the prime ' indicates the exclusion of diagonals x p = x q , p ^ q, and £(•) is an indepen¬ 
dently scattered random measure with an atomless control measure /i(-). Does l k (f n ,0 exhibits a similar 
universality phenomenon? Namely, if I k (f n ,Q converges in distribution to normal for a Gaussian £(•), does 
the convergence also hold for general class of £(•) with the same control measure /i(-)? It is known that the 
law of £(•) has to be infinitely divisible and £(•) admits the decomposition: 


m = g(b) + 


ulB(x)N(du, dx), 


(15) 


where G(-) is a Gaussian random measure on A and N(-) is an independent compensated Poisson random 
measure on R x A. See Section 5.3 of Peccati and Taqqu 0 for more details. 

One may think of adapting the approximation argument used in the proof of Theorem [2J] to the multiple 
integral case, which would involve partitioning the space A into subsets of small measure. The problem is 
that unlike the Gaussian part, the Poisson part does not scale as p(B) -j 0. To see this in the simplest 
situation, take £(-B) = P(B), where P(-) is a compensated Poisson random measure on A with control 
measure /i(-). Note that P(B) + fi(B) follows a Poisson distribution with mean p(B). Since its cumulants 
are all equal to p(B) (see (3.1.5) of Peccati and Taqqu 0 ) , and since the third moment of a centered 
random variable is equal to the third cumulant, one has E (P(B)) 3 = p.(B). This means that although we 
have the standardization 


E 


P(B)/^I(B) 


= 1, 


(16) 


we also have 


lim E 

/i(S)—^0 


P(R)/vW) 


lim E 
/^(S)—^0 


P(B) n(B )~ 3/2 > 


lim E P(B) 3 n(B)~ 3 / 2 
/*(£)-> o 


lim u(B) 1 / 2 = oo. 
o 


This will violate condition ([5]) as the partition of A becomes finer. In fact, one can show that P(B)/yf p(B) —) 
0 in probability as p(B) -j 0, which means, in view of (flTH) . that the uniform integrability of |P(£?)/ 1 Jp(B )| : 
fails. For further insights, see Rotar [12|. 
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3 Application: boundary between short and long memory 

3.1 The setting 

Bai and Taqqu jl] considered the following discrete chaos processes: 

X(n)= ^2 a(n-i 1 ,...,n-ik)e il ...ei k , (17) 

—oo<ii 

where k > 2, a(-) : Z(j_ —> R is symmetric and vanishes on the diagonals, and e^’s are i.i.d. random variables 
with mean 0 and variance 1. Note that E X(n) = 0. 

In particular, Bai and Taqqu [lj studied limit theorems for normalized partial sum process of X(n): 

x m 

Yn (*) 

where [•] means integer part, and A(N) is a suitable normalization factor. Depending on the behavior of 
a(-), the stationary process X(n) may exhibit short or long memory. 

As shown in Bai and Taqqu [1|, in the short memory case, namely when the coefficient in m satisfies 
the summability condition 


E E 

n= 1 0<ii ,...,ifc<oo 


a(h, ■ • ■, ik)a{h + n,..., i k + n) 


< oo, 


(18) 


and E|ei | 2+<5 < oo for some 8 > 0, the following central limit convergence as N —> oo holds: 

. m 

n =1 

for some a > 0, where B(t) is a standard Brownian motion. 

In the long memory case, assume that 


a(-) = g(-)L(-)l D c, 


(19) 


( 20 ) 


where 

D .— l(h? * • * Dfc) ■ 7 ^ for p / 5 } 

guarantees that a(-) vanishes on the diagonals. The function L(-) : —> R satisfied 

lim L( i) = 1, 

|i|—>-00 


( 21 ) 

( 22 ) 


and g(-) : R is the so-called generalized Hermite kernel of Class (L). 

Definition 3.1. A nonzero a.e. continuous function #(•) : —> M is called a generalized Hermite kernel of 

Class (L) (GHK(L)) if it satisfies 

1. g(-) is homogeneous with exponent a , namely, #(Ax) = \ a g(x), for all A > 0, where 


a G 



(23) 


1 In Bai and Taqqu [IJ eq. (25), L(-) is assumed to satisfy a slightly weaker condition than (122D , that is, lim^y^oo L([iVx] + 
B(AT)) = 1 for any x G and any bounded sequence B(TV) in instead of limyxii^oo T(x) = 0. Note that L([iVx] +B(iV)), 
N —>• oo, lets the argument increase in a specific band in the first quadrant, whereas T(x), ||x|| —>• oo, allows x to increase in an 
arbitrary way in the first quadrant. Here for simplicity we just assume pa . while the results stated here also hold under the 
weaker condition. 
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2. The function g(-) satisfies the bound 


| 5 (x)|< fl *(x):=c^^ 1 ...^ fc , (24) 

i=i 

with the constant c > 0 , —1 < 7 ji < — 1/2 and Yli=i 7 'jl = a: for allZ = 1 ,..., in. 


If g is a GHK(L), the following constant is well-defined (the integral is absolutely integrable) 

C g = g(xi,...,x k )g(l + xi,...,l + x k )dx 1 ...dxk, (25) 

J R* 

and C g > 0 always (Remark 3.6 of Bai and Taqqu Q). Under this setup, Theorem 6.5 of Bai and Taqqu Q 
showed that as N —> 00 , 


1 

JgH 



g(S\ X\, . . . , Sfc *^fc)l{si>a:i,...,Sfc>^fc} 


W(dxi).. .W(dxk), 


(26) 


where Wf) is the Brownian random measure, the prime ' indicates the exclusion of the diagonals x p = x q , 
p 7 ^ q, and 

H = a + ^ + 1. 

The limit in (12611 was called a generalized Hermite process which generalizes the Hermite process (see, e.g., 
Dobrushin and Major @] and Taqqu 0 ) which corresponds to the special case g(x) = x^ k ... x°^ k . 

There is, however, a boundary case which the limit theorems (fT9l) and (|26p did not cover. This boundary 
case is as follows: set as in the long memory case 


a(-) = g(-)L(-) 1 DC , (27) 

where D c is as in (12T1) . L(-) is as in (j22)l . and g is a function satisfying the assumptions in Definition 13.11 
except that instead of assuming urn the homogeneity exponent is set as a = — . 

Remark 3.2. Note that if a < — we are in the short memory regime. Indeed Proposition 5.4 of Bai 
and Taqqu @ showed that a < implies ini), and thus (fl9l) holds. So (1281) is exactly the boundary 

case between short and long memory. 


3.2 Statement of the limit theorems 


Let throughout => denote weak convergence in Skorohod space D[ 0,1] with uniform metric. We shall show 
by the criterion formulated in Theorem 12.61 that a central limit theorem holds with an extra logarithmic 
factor in the normalization: 

Theorem 3.3 (Nonlinear case). Let 


X(n)= ^2 a{n-ii,...,n-ik)e il ...e ik 

— oo<ii 

as in 0 with k > 2 and the coefficient a(-) specified as in (HP where 

k + 1 


Assume also that E|ej | 3 < 00 and C g > 0. Then 


[Nt] 


Y N (t ) := 


VNhiN ^ 

n= 1 

where a = -y/2 C g , and B(t ) is a standard Brownian motion. 


J2 X ( n ) ^ aB ( t ) 


(28) 
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Remark 3.4. Theorem 13.31 may be compared to a similar boundary case of limit theorems for nonlinear 
transform of long-memory Gaussian noise first considered in Breuer and Major [3] Theorem 1'. The proof 
there was done by a method of moments. See also Breton and Nourdin [s] who gave an alternative proof 
using the Malliavin calculus. 

Note that to apply Theorem 12.61 the process X(n) in (flTl) needs to have order k > 2. For completeness, 
we state also the corresponding result for linear process, namely, the case k = 1 in Theorem 13.31 though the 
limit theorem for linear process is classical (see,e.g., Davydov ,5§). 

Theorem 3.5 (Linear case). Let 

X(n)= ^2 a(n-i)ei, 

—oo <i<n 

where a(n) = L(n)n ~ 1 as n —> oo, and let L(n) -> c / 0, and the i.i.d. standardized noise et’s satisfy 
E|ei| 2+<5 < oo for some 6 > 0. Then as N —»• oo, 


Tiv(i) := 


1 

VN In N 


[Nt] 

71=1 


where a = \/2|c|, and B(t) is a standard Brownian motion. 


3.3 Proof of Theorem 13.31 

We first compute the asymptotic variance of the sum. 

Lemma 3.6. Let X(n) be given as in with the coefficient specified as in m and a as in \28\) . Then 
C g defined in £25 1) is non-negative. If C g > 0, then as N —> oo 


If C g =0, then 


E 


E 


N 




.71=1 


N 




2C g N\nN. 


= o(NlrtN). 


(29) 


Proof. Assume for simplicity L(-) = 1, and it is easy to extend the following arguments to the general case. 
First, since g(-) is homogeneous with exponent a = —k /2 — 1/2 by (1281) . one can write 

j(n) : = EX(n)X(0) = ^ g(ii,... ,i k )g(ii + n,... ,i k + n)l D c (ii,... ,i k ) 

0<7l,...,7fc<OO 


-k 


= n 1 ^2 g ) g l —+ 1 ,..., —+ 

n f \n n ) \n n J 

0<7i,...,7fc<OO \ / \ / 

= »-'/ g (t?A±± . + + . 

J R k + \ n n J \ n n J 

=: n~ 1 C n (g). 

Because the bounding function g* in Definition 13.11 is decreasing in every variable, the absolute of the 
integrand above is bounded by 

m 

g* (*i,.. .,x k )g* (a>i + 1,.. .,x k + 1) = c 2 ^ (x 1 + l) 73 ' 2,1 . ..a 'l il,k (x k + l) 7j2 ' fc 

3U32=1 


1 • 


.. dx k 
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which is integrable on K+ because all 7 Pi9 £ (— 1 ,— 1 / 2 ) and 


a; 7 (a; + l ) 7 dx < 00 


for any —1 < 7 , 7 ' < — 1 / 2 . 


Since g is assumed to be a.e. continuous, by the Dominated Convergence Theorem, as n —> 00 we have 


Cn(g) -+ c g 


/ g (xi,..., x k ) g (xi + 1,..., x k + 1) dxi ... dx k . 
J R fc 


Hence when C n ^ 0, one has when n > 0 

7 (n) ~ n~ 1 C g , 

and when C n = 0, one has 

7 (n) = 

We shall use the fact that if a n ~ n _1 as n —> 00 , then a « ~ InIV as N — > 00 . So when C g ^ 0, 

one has 


E 


" N 


_n=l 


AT AT—1 AT—1 

E 7(m -n 2 )=N 7(n) - E |n| 7 (n) ~ 2C g N In TV. 

wi,n 2 = l n=—AT+1 n=—AT+1 


Note that since 7 (n) ~ n~ 1 C g , the term ^ZE-iv+i I n l7( n ) ~ 2 C' g iV and is thus negligible. 

The preceding asymptotic equivalence also shows that if C g 7 ^ 0 then C g > 0 because the variance is 
non-negative. 

If C g = 0, following similar lines of argument, one gets (129[) . 

□ 


Lemma 3.7. Define the mapping (-,-)o : R 2 —> R as 


(xi,x 2 )o 


\x\ — x 2 \ ifx!^x 2 ; 

1 if x 1 = x 2 = x. 


For —1 < 71,72 < —1/2 and n\,n 2 £ {1,2,...}, we have for some constant C > 0 not depending on n\,n 2 
that 

E( n ! - P)+ ( n 2 - P)+ < C{n 1 ,n 2 )2 1+72+1 - 

pez 

Proof. For the case n\ = n 2 = n, choose C = Y2 P < n ( n — p ) 7l+72 < 00 since 71+72 < —1. When m 7 ^ n 2 , 
suppose that n\ < n 2 . Then 

00 pOO 

^(rii -p)ff{n 2 -_p)+ = ^p 7l («2 - n\ + p) 2 < / x 11 (n 2 - n\ + cc ) 72 dx 

p€% p =1 ■'° 

pOO 

= {n 2 - ni ) 7l+72+1 / y 7l (l + y) l2 dy, 

Jo 

where the integral converges. □ 

The following simple fact will be used. 

Lemma 3.8. Suppose that 7 j < —1/2 for all j = 1,..., k, k > 2, and 71 + ... + 7 *, > —k/2 — 1/2. TTien 

r 1 r 

— 2 _ 2 <7l + - — + 7r < —2 f or all r = 1 ,... ,k — 1. 

In addition, each 7 ^ > — 1, j = 1 ,.. ., k. 




Proof. The inequality 71 + ... + q r < — | is obvious. For the other inequality, suppose that 71 + ... + 7 r < 
—r/2 — 1/2 for some re {1,..., k}. Because 7 r +i,... , 7 h < —1/2, we get the contradiction: 71 + . • ■ + 7 fc < 
—r /2 - 1/2 - (k - r )/2 = -k/2 - 1 / 2 . 

Then we show by contradiction that each 7 j > —1. Suppose, e.g., 7 *, < —1. By what was just proved, 
one has 71 + ... + 7 k-i < —{k — l)/2. Thus by adding 7 *, < —1, one gets 71 + ... + 7 fc < — k/2 — 1/2, which 
contradicts the assumption. □ 


We need the following lemma, which is a consequence of Corollary 1.1 (b) of Terrin and Taqqu 141. 


Lemma 3.9. If a 1 ,..., a m , m > 2, satisfy 


m 

ai,... ,a n > — 1, ^ at + m > 1, 

i= 1 


i/ien /or any c > 0 



\xi - X 2 \ ai \X 2 - z 3 |“ 2 . . . \x m _1 - Xml 01 ™- 1 \ x m 


Xi\ am dx\ ... dx m < 00. 


(30) 


We need also the following hypercontractivity inequality for proving tightness in H[0,1] (Proposition 5.2 
of Bai and Taqqu 0 ) 

Lemma 3.10. Suppose that h £ £ 2 (Z k ) vanishing on the diagonals. Let X = J^ ieZ k h( i) ]/[p =1 k >1. If 
for some p' > p > 2 , E|ej| p < oo, then one has for some constant c pt k > 0 which does not depend on h that 


E[|AT] 1/P < Cp^-EflA'I 2 ] 1 / 2 . 


Proof of Theorem \‘3.3t Let C > 0 be a constant whose value can change from line to line. We first show 
that the finite-dimensional distributions of Y^ft) converges to those of <rB(t) using Theorem 12.61 First, the 
convergence of the covariance structure of Y^{t) to that of crB(t) follows from Lemma T3.61 the fact that for 
s < t we have 

El/v(f)b/v( S ) = \ [E Y N (t) 2 + E Y n (s) 2 - E(Y N (t) - P)v(s)) 2 ] ~ \ [El/v(t) 2 + E Y N (s) 2 - E Y N (t - s) 2 ] 

as N — > 00 , since A(n) is stationary. We now check the contraction conditions (1141) . For simplicity we set 
L(-) = 1 and t = 1. We can write 


*Zv(l)= /jv(ii,-..,u) e»! 

— oo<zi,...,ifc<+oo 


where 


/ at (* i , ■ ■ ■ ,ik) 
To simplify notation, we set 


N 


VNlnN ^ 

v n= 1 


^ ^ 9 iP' ik) ]-D c n{ii<n,...,ik<n} • 


(31) 


p = (pi,...,p r ), q= (gi,...,<?fc_ r ), 

il = (*1; • ■ • , U-r)> 12 = (ik-r+li ■ ■ ■ , *2fc-2r)j 1 = (ll j ^2)» 

and let 1 stand for a vector of l’s of suitable dimension. We also use the convention that x a = aq 1 ... x® n if 
x = (xi,..., x n ) and a = (ai,..., a n ). Let (Ex) = x\ + ... + x n if x = (aq, ..., x n ). 

Set g*{-) be as in Definition 13. II which we write by splitting x = (xi,X 2 ), where xi £ MI/ and X 2 £ R(L -r : 

m 

9* ( x i, x 2 ) = c xf J Xj 3 , f3j = (7ji,...,7jr), Vj■■ = (7j>+i,---,7jfe), (32) 

1=1 
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so that 


r k—r k 

Pi* + r h* = = c* 1 ' (33) 

2 = 1 2=1 2=1 

which we write simply as J0/3 + J0’ 7 7 = ST = a - F°r convenience, if some component Xj of x is negative, 
we set x a = 0 and hence g*(x) = 0. Then in view of (ED, ED and (l24l) . 


I(/jv *r /jv)(i)| < 


1 


N 


N\nN 


Y E 5 *( ni1 _ P’ ni1 “ 1 i)ff*(«2l - P,n 2 l - i 2 ) 


Wi,ri 2 = l P 
A/ - 71 


TVlniV 


Y Y (^i 1 - ii)^ 1 (rt 2 l ~ i 2 ) T7j ' 2 y^(nil - p)^ 1 (n 2 l - p) 


ni,n 2 =l ji,j 2 =l 

by using (l32l) . By Lemma T3.71 we have for the last sum, 


E^r 1 - p)^i(n 2 l - p)^ = Y II( ni -PnP 1,u n^ 1 ~ Pv) li2 ' v < C ( n i> 


Jo 


(£%)+(£/3 j2 )+r 


P 

Hence 


Pi,...,p r U=1 


II /n *r /jv|||fc_ 2r = Y ^ N * r / Jv )( i )] 2 


N m 


< 


c El E E (”i. > +< i »j. > + r („,i_!,)■>„ ( „ 2 i_ i 2 )i„ 


TV 2 (In TV) 2 
C 


i \m,722=1 jl,j2 = l 

N 


N 2 (h\N) 


E 


E (ni," 9 )o^ )+ ^^(rH,,»u)o^ )+( ^ )+T 


ll J2j3,i4=l rti,n2,n 3 ,ra4 = l 
^(nil - ii)^! (713! - ii) 77 ^ ^(n 2 l - h)* 1 ^ (n 4 l - i 2 ) r,:i 


< 


11 

C 


A 72 (In TV) 2 


E 


N 


Y ,n 2 )^ )+( ^ )+r (n 3 ,n 4 )r- )+( ^ )+r 


J1 >J2 J3iJ 4=1 ni,ra2,n3,ra 4 =l 


/ J + fETjj )+fc-r . (£47 ) + (£rj )+fe-r 

x (rii,n 3 ) 0 71 33 (n 2 ,n 4 )o J2 


(34) 


where we have applied again Lemma 13.71 to get the last inequality. Note that if one adds up the power 
exponents in the last expression, one gets 


(E/3*) + (£ r lj 1 ) + (E/3*) + i'^ r 1j 2 ) + (E/3*) + (E»7j a ) + (E/3*) + (Er/*) + — 4a + 2fc — —2, (35) 

by (33J) ) where the last eQU&lity of (J3^)j) is due to assumption. (j2S|) . 

Note also that by Lemma fEEl we have for r G {1,..., fe — 1} that 

- T 2~\< ( E ^)> W. W- WO < -£, 

and 

-^Y~ ~\< PVjJ, ( s »/*), ( s »?i 2 )> ( Ej 7 J4 ) < 

Let a 4 = (JO /?*) + (JO ^ja) + r be the exponent of (ni,n 2 )o in (1511) . Then 

— 1 = — r/2 — l/2^r/2 — l/2 + r<ai < — r/2 — r/2 + r = — r + r = 0. 

Define similarly a 2 ,a 3 ,a 4 for the other exponents in (1341) . which all lie strictly between —1 and 0. Hence, 
the convergence 

lim ||/jv*r fN\\l k - 2 r = °, r = 1,..., k - 1, (36) 

N —^oo 
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will follow if one shows that 

N 

sup N~ 2 "52 (ni,n 2 )o 1 (?Z2,n 3 )o 2 (n3,n4)o 3 (n4,ni)o 4 < oo, (37) 

N , 

n i ,n2,713,714=1 

where by (1551) 

— 1 < aj < 0, j = 1,... 4, op + a 2 + <23 + <24 = —2. (38) 

Let’s consider first the sum in (1371) over only distinct 7ii,..., 714 (we use the prime ' to indicate that the 
sum does not include the diagonals). In this case, 


r 

E 

l<ni ,ri2 ,ri3 ,n4<iV 


Til 

712 

Oil 

712 

713 

0.2 

713 

n 4 

03 

77-4 

?ii 

N 

N 


N 

N 


N 

N 


N 

N 


N ~ 4 


[Nxi] - [Nx 2 ] 
N 


Oil 


X 


[Nx 2 ] - [IVx 3 ] 

02 

[Nx 3 ] - [Nx 4 \ 

03 

[Nx 4 ] - [Nx 1 ] 

N 


N 


N 


/{iv 1 < Xi < i + i\r\ 


[Nxi] ^ [Nxj],Vi ^ j}dx. 


Note that for any x,y > 0, one has that |[A/rc] — [Ny]\ = n implies that \Nx — Ny\ < n + 1 which implies 
\x — y\ < {n + 1)/N, for n > 0. Then since each a < 0, we get 


sup 

N 


[Nx] - [Ny } 
N 


\x-y\ a I{[Nx] ± [iVy]} < 


sup 

| [N x] — [N y]\=n ,n€ r Lj t - 


\n) v N ) 


Hence the the sum in (1371) over distinct ni,..., 714 is bounded by 


sup 

nGZ-i- 


n -hi 
n 


— OL 


= 2 ““. 


C / \xi — x 2 \ ai \x 2 —Xz\ a2 \xz —Xi\ a3 \X 4 , —Xi\ ai dx\dx 2 dX 3 dX 4 L , 

J [ 0 , 2] 4 

which is finite due to Lemma 13.91 

Consider now the the sum in (1371) over 711 ,... ,714 with only three of them distinct. Let, for example, 
ni = 714 , and we need to show that the following 


sup TV 

N 


-2 


E K - n 2 | ai K “ n 3\ a2 \n 3 - ni|“ 3 = 

l<ni,ri2,n3<iV 


supIV 1+ai+a2+ “ 3 E 

l<ni,n2,ri3<Ar 


N 


Til 

712 

0 1 

712 

713 

02 

713 

Til 

N 

N 


N 

N 


N 

N 


TV" 3 < 00 . 


Note that (1551) entails that —2 < ay + a 2 + < —1. Then j\r 1+ “ 1+ “ 2+Q ' 3 —> 0 as N —> 00 , and the 
boundedness of the multiple sum can be established similarly as above using integral approximation and 
Lemma 13.91 

If the sum in <mzD is over 711 ,... ,714 with only two or less of them distinct, the boundedness is easily 
established through bounding all the summands by one constant, because we have the factor N~ 2 . 

So (1571) holds and thus (1551) holds, and the convergence of finite-dimensional distributions is proved. 

Now we show tightness. By Lemma [3.101 one can choose p £ (2, 3), so that by Lemma 13761 if 0 < s < t < 1, 
one has for N large enough, 


E\Y N (t) - Yiv( S )| p <C[E|l)v(t) - Y n (s) | 2 ] p/2 < C 


[Nt] - [Ns] ln([7Vf] - [IVs])l p/2 


N 


In N 


<C 


[IVf] - [Ns] 
N 


p/2-8 


where 5 > 0 is small enough so that p/ 2 — S > 1. The last inequality is true because In 2 ; is slowly varying 
as x —> 00 and so one applies the Potter’s bound (see e.g., equation (2.3.6) of Giraitis et al. 0)- Note that 
/jv(t) := [Nt]/N is a non-decreasing right continuous function on [0,1] and that Fn converges uniformly to 
F(t) := t as N —► 00 . Hence by Lemma 4.4.1 and Theorem 4.4.1 of Giraitis et al. [7|, the tightness in Z)[0,1] 
is proved. □ 
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3.4 Proof of Theorem 13.51 

Proof. Set for simplicity L(n ) = c. The covariance 7 (n) for n > 0 is 

OO OO 

7 ( 71 ) = EX(ra)X(0) = ai+nOi = c 2 ^^(* + n) - 1 i -1 . 
2=1 2=1 


Note that as n —> 00 , 

OO OO / . 

X](*+= « _1 51 (^ + 1 

i=2 i=2 ' 

The last asymptotic can be seen from: 

[ (x + l) -1 x _1 dx < f 
J 2/n J 2 


-1 


— = Tl 

n J n 


>2/n 


+ 1 


-1 /r i\ -1 

nx\ 


dx 


Inn. 


2/n 


\nx\ 


+ 1 


1 / [nx ] N 1 


nOO 

dx< {y+l)~ 1 y~ 1 dy, 
J 1/n 


where we have used the fact x—1/n< [nx]/n < x, and both the lower and upper bounds are asymptotically 
equivalent to Inn as n -> 00 . 

Hence 

7 ( 71 ) ~ c 2 n _1 Inn as n —> 00 . (39) 

So as N —> 00 , one has 

/ n \ 2 JV-l JV-l 


=N 9'( n )^ H W\l{n) 

\n—l J n=— iV+1 n=— AT+1 

]V 

~ 2c 2 fV ?r _1 Inn ~ 2c 2 N" J x _1 Inxdx ~ 2c 2 fV(lnfV) 2 . 


(40) 


n=l 


Note that by (l39l) the term X]^=-at+i I n l7( n ) = O(^V'lniV) and is thus negligible. Having obtained the 
asymptotic variance (1401) . the proof is then concluded by applying Davydov |5| Theorem 2 (though this 
theorem was stated for a linearly interpolated version of Yjv(t) in the space C[ 0, 1 ], it is straightforward to 
adapt the the proof, which consists of showing convergence of finite-dimensional distributions and establishing 
tightness by moment estimate, to establish convergence in Z?[ 0 , 1 ] with the uniform metric.) □ 

Remark 3.11. One may wonder if it is possible to get a different normalization in the nonlinear case in 
Theorem 13.31 since the normalization in the linear case in Theorem 13.51 has an extra \/ln N factor. This is 
not possible under our setting where the kernel g is homogeneous with exponent a and is bounded by a linear 
combination of products of purely power functions xj 1 ... xjf , where each 7 j < — 1/2 and 71 + ... + 7 k — cn. 

Indeed, if one wanted to get some extra logarithmic factor in the covariance 7 ( 77 ), one would set for 
example g(x±,... ,Xk) = 2 ■■■x'Jf with 7 *. = —1. But this will not achieve the stated goal. Indeed, by 
Lemma l3.8l using contradiction, we have a = 7 i + ...+ 7 fc < —kj 2 — 1/2, which falls into the short memory 
regime (see Remark l3.2D and thus the normalization is y/N as in GSD- 

Acknowledgments. This work was partially supported by the NSF grant DMS-1309009 at Boston Uni¬ 
versity. We would also like to thank the referee for his comments. 


References 

[1] S. Bai and M.S. Taqqu. Generalized Hermite processes, discrete chaos and limit theorems. Stochastic 
Processes and Their Applications, 124(4): 1710—1739, 2014. 

[2] S. Bai and M.S. Taqqu. Convergence of long-memory discrete k-th order Volterra processes. Stochastic 
Processes and Their Applications, 125(5):2026-2053, 2015. 


12 





[3] J.C. Breton and I. Nourdin. Error bounds on the non-normal approximation of hermite power variations 
of fractional brownian motion. Electronic Communications in Probability, 13:482-493, 2008. 

[4] P. Breuer and P. Major. Central limit theorems for non-linear functionals of Gaussian fields. Journal 
of Multivariate Analysis, 13(3):425-441, 1983. 

[51 Y.A. Davydov. The invariance principle for stationary processes. Theory of Probability and Its Appli¬ 
cations, 15(3):487-498, 1970. 

[6] R.L. Dobrushin and P. Major. Non-central limit theorems for non-linear functional of Gaussian fields. 
Probability Theory and Related Fields, 50(1):27 -52, 1979. 

[7] L. Giraitis, H.L. Koul, and D. Surgailis. Large Sample Inference for Long Memory Processes. World 
Scientific Publishing Company Incorporated, 2012. 

[8] I. Nourdin and G. Peccati. Normal Approximations With Malliavin Calculus: From Stein’s Method to 
Universality. Cambridge Tracts in Mathematics. Cambridge University Press, 2012. 

[9] I Nourdin, G. Peccati, and G Reinert. Invariance principles for homogeneous sums: universality of 
Gaussian Wiener chaos. The Annals of Probability, 38(5):1947-1985, 2010. 

[10] G. Peccati and M.S. Taqqu. Wiener Chaos: Moments, Cumulants and Diagrams: a Survey With 
Computer Implementation. Springer Verlag, 2011. 

[11] G. Peccati and C. Zheng. Universal Gaussian fluctuations on the discrete Poisson chaos. Bernoulli, 20 
(2):697-715, 2014. 

[12] V.I. Rotar. Limit theorems for polylinear forms. Journal of Multivariate analysis, 9(4):511-530, 1979. 

[13] M.S. Taqqu. Convergence of integrated processes of arbitrary Hermite rank. Probability Theory and 
Related Fields, 50(l):53-83, 1979. 

[14] N. Terrin and M.S. Taqqu. Power counting theorem in Euclidean space. In Random Walks, Brownian 
Motion, and Interacting Particle Systems, pages 425-440. Springer Verlag, 1991. 


13 



